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On the Scaling Law for Broadband Shock Noise Intensity in 

Supersonic Jets 


Max Kandula 1 

ASRC Aerospace , NASA Kennedy Space Center , Florida 32899, USA 


A theoretical model for the scaling of broadband shock noise intensity in supersonic jets 
was formulated on the basis of linear shock-shear wave interaction. An hypothesis has been 
postulated that the peak angle of incidence (closer to the critical angle) for the shear wave 
primarily governs the generation of sound in the interaction process rather than the noise 
generation contribution from off-peak incident angles. The proposed theory satisfactorily 
explains the well-known scaling law for the broadband shock -associated noise in supersonic 
jets. 
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Nomenclature 

= sound speed, ^ypTp 

= critical sound speed 
= specific heat at constant pressure 

= nozzle diameter 
= effective j et diameter 
= frequency 
= acoustic intensity 
= wave vector 

= wave number; also turbulent kinetic energy 

= shock-cell spacing 

= convective Mach number, u c / c ^ 

= design Mach number 
= fully expanded jet Mach number 

= transfer function for sound wave generation 
= pressure 

= reference sound pressure, 2x 1 O' 5 N/m 2 * 
= gas constant 

= jet exit Reynolds number , Pj e u je d j ! jUj e 

= sound pressure level 
= temperature 
= convective velocity 
= velocity in the streamwise direction 
= axial distance from the nozzle exit plane 
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Greek Symbols 



y = isentropic exponent 

p = dynamic viscosity 

p = density 

6 = observation angle relative to the downstream jet axis 

Subscripts 

0 = ambient fluid 

cr = critical 

e = nozzle exit 

j = jet 

ref = reference 

t = stagnation (total) 


L INTRODUCTION 

Noise from subsonic jets is mainly due to turbulent mixing, according to the theoretical model of Sir James 
Lighthill. 12 The turbulent mixing noise is essentially broadband. In perfectly expanded supersonic jets (nozzle exit 
plane pressure equals the ambient pressure), the large-scale mixing noise manifests itself primarily as Mach wave 
radiation (Ffowcs Williams 2 3 ) caused by the supersonic convection of turbulent eddies with respect to the ambient 
fluid. In imperfectly expanded supersonic jets (nozzle exit pressure different from the ambient pressure) typical of 
jet and rocket exhausts at off-design conditions, additional noise is generated (Fig. 1, see Ref. 4) in the form of 
broadband shock-associated noise (BBSN) emanating from shock-turbulence interaction (Harper-Bourne and 
Fisher 5 ) and screech tones (Powell 6 7 ) with the tonal (screech) amplitude shown to be occasioned by shock-acoustic 
wave interaction 8 . 


In imperfectly expanded supersonic jets, the rapid variation in the pressure across the nozzle exit is 
accompanied by a system of steady compression (oblique shock) and expansion waves (Fig. 2). The structure of 
these shock cells was investigated by Emden 9 , Prandtl 10 , Lord Rayleigh 11 , Pack 12 , and others. The broadband shock 
noise is of relatively high intensity, and may form a significant component of the overall jet noise, depending on the 
flow conditions. A fundamental understanding of the mechanism by which turbulence interacts with a shock wave is 
thus requisite in the analysis of the complex phenomena of shock noise generation. 

Lighthill 13 and Ribner 14,15 originally suggested that the scattering of eddies by shocks could be a strong source 
of supersonic jet noise. The importance of source coherence, however, has not been recognized, so that only 
incoherent and randomly scattered sound waves had been predicted without the peak frequency and directivity 
relationships. It was Harper-Bourne and Fisher 5 who first identified the detailed characteristics of BBSN with the aid 
of measurements from conical (convergent) nozzles, and indicated the importance of source coherence. The 
characteristics of shock noise were also reviewed and discussed in Fisher et al. 16 Howe and Ffowcs Williams 17 also 
considered that the primary source of broadband shock-associated noise is a consequence of the interaction between 
large scale structures (turbulence) and the shock structure. 

Computing shock noise intensity in supersonic jets from first principles (on the basis of shock-turbulence 
interaction) is very difficult. The nature of the relevant noise sources is not well understood (Krothapalli 18 ). This 
situation is exemplified by the fact that the theories of both Lighthill 13 and of Ribner 19,20 produce shock noise 
intensity scaling considerably different from that indicated by the measurements. 
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It is the purpose of this work to investigate the scaling of shock noise intensity from fundamental 
considerations. It is demonstrated here that the scattering of turbulence by the leading shock wave is related to the 
measured shock noise intensity scaling. Flight effects are excluded from consideration here. Also screech effects are 
not relevant to this investigation. 


II MEASUREMENTS AND CHARACTERISTICS OF BROADBAND SHOCK NOISE 


A. Intensity of Shock Noise 

Harper-Bourne and Fisher 5 were the first to identify significant features of shock-noise in considerable detail 
based on their static jet measurements from conical nozzles. The intensity of BBSN is shown to be primarily a 
function of the nozzle (operating) pressure ratio (NPR= p t / p 0 ). For a given radiation direction, the measured 
overall sound intensity / has been observed to scale as (Fig. 3) 

Io c/? 4 (la) 

where = [m^ -\ j /2 (lb) 


with the isentropic relation between p t / /? 0 and M y expressed by 
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In the preceding relations, the quantity My represents the fully expanded jet Mach number, and the parameter 

P characterizes the pressure jump across a normal shock at approach Mach number My . Harper-Bourne and Fisher 
(1973) found that the parameter p correlated BBSN quite well up certain values of P (or NPR), say 0.5 < P < 1.2 . 
At large NPR or /? , the data begin to deviate from this law because of the presence of a Mach disc, which 
significantly alters the shock-cell structure. As the Mach disc forms, the large central portion of subsonic flow 
formed downstream of the Mach disc considerably reduces the noise generation. The data also reveal that at high 
P the turbulent mixing noise level is much lower than the underexpanded noise levels. As p decreases, the mixing 
noise contribution relative to the total noise becomes increasingly significant. 

Experiments by Tanna 21 and of Seiner and Norum 22 provided further insight into the characteristics of shock 
noise. These data include measurements form convergent-divergent (C-D) nozzles, and covered a broad range of jet 
conditions (NPR and jet temperature ratio). Both Tanna’s data 21 , covering p < 1 (m y <1.41, or p R l /? 0 < 3.5), and 

the data of Seiner and Norum 22 (covering M d = 1.5, and 2 and My = 1 to 2.37 or p= 0 to 2.15) suggest trends 
similar to those indicated by the data of Harper-Boumer and Fisher 5 to the extent that the overall intensity of shock- 
associated noise is principally a function of jet pressure ratio, scales as / oc /? 4 , and is independent of jet temperature 
ratio (efflux temperature) and emission angle. 

The data by Krothapalli et al. 18 for broadband shock noise for My in the range of 1.24 to 1.66 suggest that the 

shock noise intensity follows the P 4 dependence for both the stationary ambient ( M^ = 0 ) and in forward flight 
(M^ = 0.32 ). 
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B. Spectra and Directivity of Shock Noise 

Directivity 

Measurements by Tanna 21 indicate that for all values of Tj I 7q considered, the directivities from underexpanded 

jets in the forward arc are essentially flat, indicating that the sound radiated by the presence of shocks in a jet flow is 
fairly omnidirectional (Fig 4). 

Detailed acoustic measurements by Norum and Seiner 23 suggest that the shock noise intensity is not necessarily 
omnidirectional. The shock noise is fairly directional at lower fi values and approaches omni directionality only at 
large values of P . Both the source directivity and the Doppler shift effect may contribute to the far-field directivity 
(Pao and Seiner 24 ). 

According to Krothapalli et al. 18 , in cold jets, the shock-associated noise dominates over the turbulent mixing 
noise except in the aft quadrant. At elevated temperature, the turbulent mixing noise increases, and the shock- 
associated noise remains unaltered. 

Spectra 

Fig. 5 displays the spectral characteristics of BBSN according to the data of Tanna 21 , for various angles of 
observations. The variation of peak frequency (which represents an important characteristic) with the angle of 
observation is demonstrated. Pao and Seiner 24 indicate that the power spectral density (dB/Hz) increases as 

&> 4 below the peak frequency, and decays as co~ 2 beyond the peak frequency. 

According to the measurements by Krothapalli et al. 18 for broadband shock noise for Mj from 1.24 to 1.66, the 

broadband spectral peak is dominant in the region of 50 < 0 < 1 10 deg. Unlike the screech tone, its peak frequency 
is a function of the direction of radiation with the lower frequency being near the jet inlet direction ( 0 = 0 ). 


Ill MODELS FOR BROADBAND SHOCK NOISE 

The semi-empirical model of Harper Bourne and Fisher 5 forms an important contribution to the theory of BBSN. 
Alternative models were later proposed by Howe and Ffowcs Williams 17 and Tam and coworkers. 25-27 Howe and 
Ffowcs-Williams 17 proposed a theoretical model of the BBSN based on a multiple scattering approach. Tam and 
coworkers 26,27 proposed a stochastic model theory for BBSN that showed good agreement with experimental data. 
These models deal with shock noise intensity, spectra and directionality. 

A. Model of Harper-Bourne and Fisher 

Harper-Bourne and Fisher 5 proposed a semi-empirical physical model as an extension of Powell’s model for 
screech. 6,7 Central to their investigation is the intensity scaling provided by Eq. (la), namely 

7oc/? 4 

.They approximated (idealized) shock-noise sources as a series of point sources located at shock-turbulence 
interaction locations. The scaling of the spectral peak frequency with shock-cell length, and the correlated nature of 
shock noise emission from consecutive shock cells were investigated. It was shown (demonstrated) that if the 
turbulence maintained its coherence as it convected through multiple shock cells, positive (constructive) interference 
from the sources would result in strong radiation in the upstream direction. 

The predicted (and observed) peak frequency was expressed by (see also Fisher et al. 16 ) 
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(2a) 


f -£ 

Jp L(\ - M c cos 0) ’ 

where w c is the appropriate convection velocity, L the shock-cell spacing, M c (=u c /c OQ ) the convective Mach 
number, and 6 the observation angle (direction of noise radiation) relative to the downstream-directed jet axis (jet 
flow direction). The shock-cell separation L is expressed by 

L = c\d jfl (2b) 

where dj is the nozzle exit diameter, and q is a constant of proportionality. This relation was theoretically derived 

by Pack 12 . From the measurements, an average value for the constant for the shock-cell array (about eight shock 
waves) is found to be 1 . 1 . No harmonics are observed in the noise spectra, 

A general prediction method for the sound spectra was also developed, requiring input concerning the source 
spectral density and correlation coefficients, which can be obtained from experiments. The predictions from the 
model are shown to satisfactorily agree with the data. . . 

B. Model of Howe and Ffowcs-Williams 

Howe and Ffowcs-Williams 17 proposed a theoretical model for the BBSN, considering the interaction between 
turbulent velocity fluctuations and the shock-cell system of an imperfectly expanded supersonic jet. Their procedure 
constitutes a generalization (application) of Lighthill’s acoustic analogy 1 of aerodynamic sound to the scattering of a 
sound wave by turbulence. According to this model, energy is extracted (abstracted) from the ordered nonuniform 
mean flow of the jet and redistributed into random scattered disturbances. The peak of the noise spectrum is 
associated with the coherent scattering of sound by the array of shock cells, with additional sound waves produced 
through multiple scattering forming a wide frequency range of the broadband noise. Isothermal jet is considered, and 
formulas were derived for the spectra and directivity of BBSN. 

With regard to the noise intensity, the model reveals that the predicted sound intensity varies as 

/ *fi 5 6 (3) 

This relation is at variance (deviates from) the experimentally observed law (Eq. la)., showing a fourth power 
dependence on p It is pointed out that increasing/decreasing the shear layer width tends to decrease/increase the 
effective exponent of the predicted power law dependence on p (Howe and Ffowcs Williams 17 ). 

The predicted sound spectrum is characterized by a sequence of peaks produced on account of coherent 
scattering from successive shock cell. The peak frequency is theoretically found to be consistent with that expressed 
by Eq. (2), except that the constant is shown to be 1.31. The smaller peaks of the predicted multiple-peaked 
spectrum are not harmonically related. The directivity of the broadband shock noise is somewhat (essentially) 
uniform over a wide range of angles in the forward arc. 

C. Model of Tam and Coworkers 

Tam and Tanna 25 presented a shock noise model which considers (weak but coherent) interaction between the 
downstream propagating instability waves and periodic shock cell system (turbulent eddies regarded as instability 
waves or large scale turbulent flow structures in the jet mixing layer). Generalizing Harper Bourne and Fisher’s 

p 4 intensity scaling to convergent-divergent nozzles, they proposed a correlation of the form 

I X ( M j-M 2 d f (4) 
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where M d refers to the nozzle design Mach number. Eq. (4) thus applies to both convergent (M l d =1) and 
convergent-divergent nozzles. This formula is valid only for weakly imperfectly expanded jets. For convergent 
nozzles, it reduces to scaling relation of Harper-Bourne and Fisher 5 . Data for a C-D nozzle for 
M d = 1 .67 ( M j =1.1 to 2.0) were also presented, including frequency spectra 

A frequency-directivity relation, equivalent to Eq. (2), was derived., The peak frequency is shown to be 
expressed by 


/ = — 7 — =* t (5a) 

Lj (l - M c cos 9) 

where v Lj = — = jt{Mj - if d e / //j (5b) 

l i 

with ju\ =2.40483, and d e is the effective jet diameter. The quantity d e !dj is a function of M d ,Mj and yby 

virtue of conservation of mass flux and isentropic relations (Tam and Tanna 25 ). The above expression for the peak 
frequency is similar to that of Harper-Bourne and Fisher 5 , as indicated by Eq. (2). With this model, it was suggested 
that shock-noise should radiate strongly at upstream angles as Mach waves. 

Building on the work of Tam and Tanna 25 , Tam 26 formulated a stochastic model theory of the BBSN of 
axisymmetric supersonic jets by considering the dynamics of weakly nonlinear interaction between linear instability 
waves (superposition of intrinsic or instability wave modes of the flow describing the spatially and temporally 
evolving large turbulent structures in the mixing layer) and shock-cell structures (with the standing wave pattern 
inside the wave guide enclosed by the plume boundary, obtained by multiple-scales expansion). In view of the 
solution complexity, a final solution was not arrived at; instead a semi-empirical (less general) shock noise model, 
valid for slightly imperfectly expanded supersonic jets. Formulas for the nearfield and the farfield spectra were 
provided by aid of physical reasoning, scaling arguments and data correlation. The results satisfactorily explain in 
many important aspects of the noise spectra and its directivity pattern. Typically an increase in the spectral peak 
associated with the BBSN is attributed to the convective amplification of the sources. 

The stochastic model theory of BBSN by Tam 26 was extended by Tam 27 to moderately imperfectly expanded 
jets with the aid of empirical modifications to the amplitude of the waveguide modes of the shock cell. A 
mathematical theory for the interaction between large turbulent structures and quasi-periodic shock cells is thereby 
proposed. 


IV FUNDAMENTAL STUDIES ON SHOCK-TURBULENCE INTERACTION 
A. Theoretical Models 
1 . Linear Theories 


Broadly speaking, the decomposition of a general fluctuation into acoustic, vorticity and entropy waves is well- 
known (Kovasznay 28 ). In general any plane wave (acoustic, vorticity/shear, or entropy) interacting with a shock 
undergoes transformation, and at the same time generates the other two waves (Zang et al. 29 ). For example, Fig. 6 
shows a schematic of an interaction of a shock wave with an incident vorticity wave. Linear analyses of a single 
wave (shear/vorticity, acoustic, or entropy) interaction with a shock wave were carried out by Blokhintzev 30 , 
Burger 31 , Ribner 14,15,19,20 , Moore 32 , McKenzie and Westphal 33 . 

According to the linear theory, for sufficiently high angles of incidence for the wave ahead of the shock, the 
incident wave vector k defined by 
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= 2 k! X = col c 


( 6 ) 


t.|K| 

(where X is the acoustic wavelength, cthe sound speed, and co the circular frequency) has a nonzero imaginary part. 
Under such circumstances the refracted (or generated) acoustic wave is not an infinite plane wave; instead, it 
exhibits an exponential decay as it propagates downstream behind the shock. The incidence angle that separates the 
plane wave acoustic response from the decaying ones is termed the critical angle . The critical angle is close to 90 
deg. for incident acoustic waves, and roughly 60 deg. for incident vorticity and entropy waves (Zang et al. 29 ). Linear 
theory predicts that most transmission and generation coefficients are peaked near the critical angle (evanescent 
acoustic response?). From a theoretical point of view, the actual transmission/generation coefficients are 
independent of the incident wavelength in the linear limit (Zang et al. 29 ). 

With regard to broadband shock noise, we are primarily concerned here with shock-turbulence interactions. A 
turbulent velocity field can be represented as a superposition or spectrum of elementary waves distributed among all 
orientations and wavelengths in accordance with Fourier’s integral theorem. The waves are transverse for weak 
turbulence because of the constraint of incompressibility (even though convected at high speed). Thus a single wave 
can be interpreted physically as a plane sinusoidal wave of shearing motion (Batchelor 34 ). According to linear 
Interaction Analysis (LIA), the vorticity waves incident at angles beyond a critical angle 


e c =0 c( M l) (7) 

generate acoustic waves which decay as they propagate downstream. In eq. (7), Mi refers to Mach number upstream 
of the shock 

Lighthill 3 and Ribner 14,15,19 conducted theoretical analysis on acoustic noise generation by shock 
wave/turbulence interaction. In both Ribner’ s and Lighthill ’s theories the turbulence is treated in effect as a frozen 
spatial pattern with neglect of temporal fluctuations. 

Ribner ’s Analysis 

Ribner 13 studied in detail the interaction between a vorticity wave and a shock wave. Ribner 14,18 extended this 
analysis to consider a spectrum of incident vorticity waves (in three dimensions) and computed, for an isotropic 
incident spectrum, detailed statistics of the downstream flowfield with emphasis on the generated noise. The basic 
building blocks of Ribner’s linear theory are oblique plane sinusoidal waves of vorticity (shear waves), see Fig. 6. 
These represent single spectral (monochromatic) composed of (in 3-D) an instantaneous snapshot of arbitrary flow. 
The waves are considered to interact independently with the shock, and then the waves are superposed to represent 
turbulence upstream and downstream of the shock. The detailed statistical formalism was worked out in Ribner 15 , 
and partly summarized by Ribner 20 . 

The mean spectral sound pressure is expressed by 15,19 

P " 2 = T| P(&)\ 2 [uu]d 3 k (8) 

—00 

where P(&) is the transfer function for sound wave generation, and [«w]the longitudinal spectral intensity. The 

initial turbulence is now restricted to be isotropic, so that its longitudinal spectral density has the general form 
(Batchelor 34 ) 


[uu\ = k 2 F(k) cos 2 0 

in spherical coordinates in wave number space, where F{k) is an arbitrary function of k . 


(9) 


For an isotropic turbulence, the acoustic energy flux (per unit area normal to the shock) is finally expressed by 
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(10a) 


*22 ^ /2 

/ ac = - £- 2 -? . f |P((9)| 2 cos 3 <9(l + M 2 sin 6 »")(a/ 2 + sin 9")d0 

2 pi c i i 

°cr 

2 

or — — = ^ f |/ > (^)) 2 cos 3 ^(l H- A^ 2 sin <9"X^2 +sin&”)d@ (10b) 

7 turb 5 pip 2 c* c 2 U\ ^ 

where the turbulence intensity 7 turb is defined by 

/turb=|^lt/l^(«l 2 ) ' (10c) 

Here 6" denotes the angle of the refracted shear wave, U the flow speed upstream of the shock, c* the critical sound 
speed, and the ambient pressure. It is necessary to obtain the transfer function P(o) and the incident shear wave 
inclination 6" - f(0) . These quantities are rather cumbersome functions of other functions obtained by Ribner 14 
and evaluated and tabulated in Ribner 15 along with the integration limit O cr . Calculations of the linear theory were 
performed for an upstream Mach number range of 1 < M\ < 10 . 

The sound pressure level (SPL) is given by 20 

SPL = 20\og(j7/ Pref ) (11) 

where /? ref is the reference sound pressure (2xlO" 5 N/m 2 ). For one percent turbulence, the post shock noise level is 

predicted to exceed 140 dB for all preshock Mach numbers above 1.05. The velocity components of the post-shock 
turbulence would be amplified as much as 45 % relative to preshock levels. 

Lighthills ’s Theory 

Lighthill 13 considered the generation of sound due to the interaction of turbulence with very weak shock waves 
(acoustic-like waves), by aid of his general theory of sound generated aerodynamically (Lighthill 1,2 ). The weak 
shock is represented by an acoustic step function.. In Lighthilfs theory the assumptions are more restrictive than in 
Ribner’ s analysis in the sense that both the shock and the turbulence are weak. As a result, the rippling motion of the 
shock as well as the differences in the turbulence intensity across the shock are suppressed. The ratio of freely 
scattered acoustic energy to the kinetic energy of turbulence traversed by the shock wave is expressed relative to a 
frame moving with the fluid , whereas Ribner’s analysis deals with a frame attached to the shock. 

For a direct comparison, Ribner 19 converted the results of Lighthill 13 to the shock-fixed reference frame as 
follows: 


7^4 

1 turb 


2175 -6<z> + ^ 3 L 2 


(12a) 


where 


' 3, ^ 


10 + 6 * 


£ = 


'a-.' 

Pi 




(12b) 


with the density ratio across the shock expressed by the Rankine-Hugoniot relation for a normal shock. 
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Fig. 7 shows a comparison of the scattered sound intensity (in SPL) between Ribner’s result (Ribner 19 ) and that 
of Lighthill 13 , as presented by Ribner 19 . Significant discrepancy is noted between the two results. A critical' 
discussion of this comparison is provided by Ribner 19 . (mention the scale for Ribner’s and Lighthill’s results) 

2. Nonlinear Euler Simulations 


Since the shock weakens as the Mach number tends to unity, the shock front will undergo greater distortions 
from an incident wave of fixed amplitude. Thus, nonlinear effects ought to be increasingly important for lower Mach 
numbers (Zhang et al. 29 ). Zhang et al. 29 validated the linear analysis of McKenzie and Westphal 33 by comparisons 
with their numerical solution of nonlinear 2D Euler equations. Although restricted in terms of the incident angle of 
the disturbance, it was shown that the linear analysis was valid over a surprisingly large range of shock strengths and 
disturbance amplitudes (Fig. 8). The comparisons suggest that the linear theory is fairly accurate for a wide range of 
incident angles up to the critical angle. Although comparisons were made for both the incident and the vorticity 
waves, only the comparisons for the incident vorticity waves are indicated in Fig. 8. 

3. DNS Simulations 


The simplest circumstance in which turbulence interacts with a shock wave is the case of isotropic turbulence 
interacting with a normal shock (transverse vorticity amplification). Lee et al. 35 performed DNS simulation of the 
interaction of 3-D isotropic turbulence (0.0 6<M t <0.11) where M t is the fluctuation Mach number (i.e., rms 

velocity/mean speed of sound upstream of the shock wave) with a normal shock ( 1 .05 < M x <1.2) where M x is the 
mean upstream Mach number. Detailed comparisons of DNS results to Ribner’s linear analysis 15,20 were made. 
Amplification (of turbulence) is more pronounced at the large wave numbers, which is consistent with the linear 
analyses. The turbulent kinetic energy is amplified across the shock wave, but this amplification tends to saturate 
beyond M - 3 . DNS results for the amplification of kinetic energy of turbulence across a normal shock wave vs. 
linear theory is good. The rapid variation of k immediately behind the shock is connected with the decay of the 
acoustic energy in this region. The unsteady motion of the shock wave is found to reduce the amplification of the 
turbulent kinetic energy. 

Subsequently stronger shock strengths (M x - 2, 3; M t =0Al) were studied (Lee et al. 1997). DNS calculations 
by Lee et al. 36 agree with Ribner’s theory for the measure of anisotropy downstream of the shock as u 2 / v 2 =0.87 at 
M=3 in agreement with Ribner’s theory. The DNS studies also suggest that the isotropic turbulence becomes 
axisymmetric downstream of the shock. Numerical simulations by Rotman 37 and DNS calculations by Lee et al. 34,35 
show that the vorticity amplification predictions are in good agreement with the linear theory. 

Fig. 9 shows a comparison of amplification of turbulent kinetic energy ( k 2 /k x ) in shock-turbulence 
interaction, as predicted by linear theory and DNS (Mahesh et al. 38 , Lee et al. 36 ). This figure is adapted from Sinha et 
al. 39 ). Excellent agreement is noticed between the two computations up to M x = 3 , beyond which the linear theory 
suggests an asymptotic trend in the amplification. 

Although DNS solutions provide the most accurate representation of the shock/turbulence interaction, they 
seem to be impractical for conditions involving strong shock waves and high Reynolds number turbulence on 
account of resolution requirements of shock waves and turbulence. . 

B. Experimental Data 

Experimental data suggest that in general compression enhances turbulence and expansion suppresses it. 
Measurements by Barre et al. 40 at M\ = 3 suggest that the shock wave increases the longitudinal fluctuating velocity 
in agreement with Ribner’s theory (Ribner 19 ). As indicated by Ribner 41 , the measured amplification ratio of mean 
square longitudinal component of turbulence velocity ( u 2 / u{ ) is close to the theoretical value of about 1.5 as 
predicted by Ribner’s theory 14,15 at M\ = 3 (Fig. 10a).. 
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With regard to the anisotropy of turbulence downstream of the shock, the DNS calculations of Lee et al. 34,35 and 
the linear theory of Ribner 15 suggest that the anisotropy downstream of the shock u 2^2 = 0.87 at A/] =3 . This 
result considerably deviates from the measured value of 1.5 (Barre et al. 40 ). as seen in Fig. 10b. Linear theory 
predicts a value of u 2 !u{ in excess of unity only for upstream Mach numbers below 2, with the highest value shown 
to be about 1.1 (at M\ =1.25 ). As indicated by Ribner 41 , this discrepancy between the theory and the experiment 
remains puzzling. 


IV PROPOSED THEORY 

The discrepancy between the theories of Lighthill 13 and of Ribner 14 in comparison with the experimental data 
for the scaling of shock noise intensity (as evident from Fig. 7) requires further investigation. This discrepancy is 
attributed to the fact that in their theories the turbulence is treated effectively as a frozen spatial pattern without 
regard to the temporal fluctuations. There is thus a deficiency in applying linear theory to real turbulence, which 
consists of transient phenomena and not steady plane waves (Zhang et al. 29 ). Also, three-dimensional simulation is 
needed to accommodate vortex stretching (Zhang et al. 29 ). 

The irregularity and disorderliness characterizing turbulence involve the impermanance of the various 
frequencies and of the various periodicities and scale (Hinze 42 ). In view of these circumstances these circumstances 
it is plausible that the peak angle of incidence is representative of the shock-shear wave interaction insofar as the 
scaling of the BBSN is concerned. Accordingly it is postulated here that the generation of sound at the peak 
incidence governs the generation of sound. Also it is assumed that the interaction of turbulence with the leading 
shock cell forms the maximum contribution to the intensity sound, and that the sound contribution due to the 
interactions at the subsequent shock cells is of secondary nature (subsidiary importance). 

With the above postulate, the linear acoustic response (acoustic pressure rise) for shock-vortex interaction 
(vorticity waves incident on a shock) is shown in Fig. 12 for various upstream Mach numbers, as given by Ribner’ s 
theory (Ribner 14,15 ). The results point out that the peak angle of incidence and the associated acoustic response varies 
with the Mach number. 


IV. RESULTS AND DISCUSSION 

Based on the foregoing premise, the intensity of BBSN taken at the peak incidence angle is plotted as function 
of p (Fig. 10). This result reveals that the scaling of intensity very nearly varies as p 4 for a range of p between 0.2 
and 2.0. In this range the present theory yields 

/OC /? 42 (13) 

Beyond this range there is seen a change in slope in the intensity variation. 

A direct comparison of the scaling based on the proposed theory and the experimental data of Tanna 20 is 
presented in Fig. 13. The theory substantially agrees with the data in the range of 0.3 < p < 1.0 . It is known that 
beyond about /?='l, a Mach disc is formed, which alters the shock-cell structure. The large central portion of 
subsonic flow that develops downstream of the Mach disc considerably diminishes the noise generation. 

The satisfactory explanation of the /? 4 scaling law by the proposed theory suggests that the hypothesis of peak 

incidence angle for the generation of sound by shock-vorticity interaction is plausible. This forms an important 
contribution of the present work. 

The determination of the directionality effects and spectral distribution of the BBSN are outside the scope of the 
present investigation, which is mainly concerned with the scaling law for broadband shock noise intensity. The fact 


10 

American Institute of Aeronautics and Astronautics 



that only a single shock-cell/vortex interaction is considered here indicates that the shock noise intensity obtained by 
the present formulation is essentially ominidirectional. 

V. CONCLUSION 

A theory is proposed for the scaling of the broadband shock-associated noise in supersonic jets considering 
linear interaction between the shock wave and the vorticity wave on the basis of peak incidence angle for the 
turbulence. The hypothesis that the generation of sound at peak incidence angle is important is shown to 
satisfactorily describe the experimental scaling law for the broadband shock-associated noise intensity. 
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Fig. 1 A typical narrowband farfield shock noise spectrum (adapted from Seiner 4 ). 
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Fig. 2 Shock cell structure in imperfectly expanded supersonic jets. 
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Fig. 3 Intensity of broadband noise at 90 deg. to jet axis (from Fisher et al. 16 ) 
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Fig. 4a Directivity of broadband shock noise (from Tanna 21 ) at /? = 0.94 
and Tj i T 0 =0.73. 
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Fig. 4b Directivity of broadband shock noise (from Tanna 21 ) at /? = 0.94 
and Tj/T 0 =2.272 . 
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Fig. 5 Spectral character of broadband shock noise (from Tanna 21 ). 



Fig. 6 Interaction of a shear wave with shock wave (Ribner 19 ). 
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Fig. 7 Intensity of broadband noise according to the theories of Ribner 19 and Lighthill 13 . 



Fig. 8 Dependence of acoustic response to vorticity waves incident on a Mach 8 shock (solid lines: linear theory, 
circles from nonlinear Euler simulations; from Zhang et al 29 ). 
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Fig. 9 Amplification of tumbulence kinetic energy in shock-turbulence interaction, 
as obtained by linear theory and DNS (Mahesh et al. 37 ) 
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Fig. 10a Measured amplification of longitudinal velocity fluctuation across a shock 
in shock-turbulence interaction (Barre et al. 40 ). 
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Fig. 12 Intensity of broadband noise according to the theory of maximum incidence angle 



Fig. 13 Comparison of the present theory and the data of Tanna 21 for the intensity of the 
broadband shock noise. 
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